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Abstract 



We investigate the influence of tlie pfiase-mismatcli in tfie emission 

rate of photon-pairs and in their degree of entanglement in nondegen- 

^ I erate parametric processes. The Hamiltonean describing this processes 

^O ' contains explicitly the phase-mismatch and is solved exactly. The time 

dependent bosonic operators obey to generalized Boguolibov transfor- 
f^ i mation. We show that the average number of photons produced in 

^ I the interaction, the degree of squeezing as well as the degree of entan- 

^D ■ glement decrease with increasing phase-mismatch. 



PACS number(s) 42.65.Lm, 03.67.a, 42.50-p. 



1 Introduction 

A significant number of Quantum Computation and Quantum Information 
protocols use the quadrature amplitudes of the quantized electromagnetic 
field [1], [2]. The essence of all these protocols rely in the possibility of 
producing entangled photons in optical nonlinear mediums [3]. 

In particular, the nondegenerate parametric amplifier (NPA) with in- 
jected signal, and the spontaneous parametric downconversion (SPD) gener- 
ate squeezed states of light [H |5l [6] and are paradigmatic processes in the 
realm of the continuous variable (CV) entanglement [7]. 

The efficiency of any nonlinear optical process depends strongly on this 
wave-vector mis- match. The emission is brightest if the various fields are 
coherent over the full length of the 8 

Dispersion in the nonlinear material between pump, signal, and idler 
waves hereby plays an inhibiting role and should be eliminated. This can be 
achieved with the use of birefringent nonlinear crystals. In such crystals per- 
fect phase matching can be attained by appropriately orienting the crystal 
axes, and the wave vectors and polarization vectors of the input fields. How- 
ever, it is often the case that the strongest nonlinearities cannot be addressed 
in this man- ner. Consequently, considerable effort has been devoted to de- 
velop materials that combine perfect phase matching with a strong nonlinear 
response, resulting in a multitude of quasi- phase-matched and periodically- 
poled materials [88, 89]. 

It is then worth to inquire how the quality of entanglement achieved in a 
parametric amplifier and in a spontaneous parametric downconversion can be 
affected when the phase matching condition cannot be completely achieved 
in actual experiments. It is known from classical nonlinear optics [S] that 
the existence of a phase mismatch (i.e. a nonzero Ak), that occurs during 
the propagation in the nonlinear medium, modify the power of the signal 
and of the idler. Therefore, absence of perfect phase matching will certainly 
affect efficient amplification, degree of squeezing and entanglement of the 
interacting modes. 

In this paper we generalize the usual analysis of the nondegenerate para- 
metric processes [lOl [9] to take into account the infiuence of the phase mis- 
match in the degree of entanglement. In the classical version of this problem 
[H] the phase-mismatch is mathematically equivalent to a detuning in the 
frequency of the interacting modes. The corresponding quantized version is 
described by an extension of the usual Hamiltonean [3 [ID] which is no more 



time independent in the interaction picture. We compute the exact solution 
of the Heisenberg equations of motion to find the time dependent bosonic 
operators that exphcitly depend on the phase-mismatch. 

This model also describes nondegenerate parametric interactions of waves 
with finite spectral bandwidth |TT] ,this is, the energy conservation rule is 
obeyed within the imprecision of the finite bandwidth of the pump field. 

The characterization of the entanglement will be done using the well 
stablished techniques of continuous variables Gaussian states developed for 
bipartite systems [12], this is, we use the formalism of the covariance matrix 
to compute the Von Neumann entropy of each mode as well as the logarithmic 
negativity of the composite system, for different initial quantum states. In 
all cases we prove that the degree of entanglement decreases with increasing 
mismatch or with finite bandwidth. 

The plan of the paper is the following: In Sec. II we introduce the Hamil- 
tonean of the system with phase-mismatch, solve the Heisenberg equations of 
motion for the bosonic operators and compute the average number of gener- 
ated photons. In Sec. Ill we compute the covariance matrix of the system for 
different initial states and determine in each case the degree of entanglement 
of the system. Finally in Sec. IV we present some conclusions. 

2 System dynamics 

In the NPA a pump mode with frequency 2uj and wavevector 2k interacts 
in a nonlinear medium with a mode of frequency Ui and wavevector ki (the 
signal ) which grows. In the process an auxiliary mode of frequency uj2 and 
wavevector /c2 ( the idler ) is created. The SPDC differs from the NPA 
because the only input to the nonlinear medium is the pump mode and 
downconversion to the lower-frequency modes Ui and U2 is spontaneous. We 
consider that the incident pump mode is intense and can be treated classically 
as a field of complex amplitude A = |^|e~^*'^*, only modes 1 and 2 are 
described by the bosonic operators cii and 0,2- 

In both cases the three modes obey the energy conservation rule 2u = 
Ui + U2 and the phase mismatch Ak = 2k — ki — k2 y^ 0. The Hamiltonean 
describing both processes and that generalizes [H [10] , is given by [8] 

H = huid\ai + hu2ald2 - ng{e-^'^'^+^'>'d\dl + e2*(^+'5)*aia2) (1) 

where the interaction time t is the propagation time in the nonlinear medium 




therefore proportional to the travelled distance z. The phase factor e^*''* 
with 5 > results from the mismatch along z and the coupling constant g 
is proportional to the second order susceptibility of the medium and to the 
amplitude of the pump. 

The Hamiltonean ([1]) can also describe the interaction between a classical 
wave of central frequency 2a; and bandwidth 25, with two quantized modes 
of frequencies cji and lo^^ and such that the energy conservation rule 2a; = 
uj\ + uji is obeyed. 

The Heisenberg equations of motion for modes 1 and 2 are 

jointly with their hermitian conjugates, which are decoupled from these two. 
When the mismatch is zero the integration of these equations is done in the 
interaction picture. However, this usual approach doesn't work in the case of 
a nonzero mismatch because it doesn't remove the time factor e~^*^* from the 
Heisenberg equations. In order to obtain an autonomous system of equations 
we are going to define new annhilation operators hx and h^^ by the following 
canonical transformation 

(3) 

Physically this corresponds to move to a reference frame that rotates with 
the frequency (w + 5). The Heisenberg equations for the new operators h\ 
and 62 cire then given by 






—ig i{u 




(4) 



They are linear time independent differential equations and can be readily 
integrated. The only solution of eq. (jl]) with physical interest corresponds to 
the condition 6 < g where the mismatch is small compared with the coupling 
parameter. By inverting eq.Q and using the solution to eq.(|l]) we obtain 
finally the time dependent bosonic operators 

(5) 



al{T) = e^™2(s/) 



W 



;S{xT)aio + C{xt) - —-prS{xT 



(6) 
where t = gt is an effective dimensionless interaction time, y = d/g is the 
dimensionless mismatch parameter with values in the interval < y < 1, the 
effective frequency of oscillation of mode j is u}j{y) = ojj/g + y. 

The system of equations ([S]) and ([H]) is a kind of generalized Boguolihov 
transformation where C{xt) = cosh.{xgt), S{xt) = sinh(x5't) are the hyper- 
bohc functions, x = (1 — y^Y^"^ and aio = ai(0) and 020 = 02(0) ^ire the 
bosonic operators at time t = 0. 

Equations ([5]) and ([6]) extend the well known result ^ to take into account 
the influence of the phase-mismatch. The corresponding bosonic operators 
derived in [9] are obtained by making y = in the last equations. 

The argument of the hyperbolic functions is the squeezing parameter, 

r{y)=gt{\-y')'l\ (7) 

which is a generalization of the usual squeezing parameter r = gt. It de- 
pends not only on the interaction time t, but also on the mismatch. For 
a given interaction time the intensity of the squeezing parameter decreases 
with increasing mismatch. 

From the Heisenberg equations ([2]) we obtain the following integral of 
motion 

(rii(t)) - (ri2(t)) = fiio -^20 (8) 

which is a generalization, for nonzero mismatch, of the analog integral of 
motion derived in [9]. We have introduced in the previous equation the initial 
mean number of photons nio and n2o- This integral of motion says that each 
time the pump creates a photon in mode (1) it also creates a photon in mode 
(2) therefore, the difference of photons between these modes at any time t is 
equal to their difference at the initial time t = and that it is independent 
of the phase-mismatch. For uncorrelated input vacuum states the average 
number of photons in each mode is initially, tt-iq = ?t-2o = 0, jointly with 
equality ([S]) it determines that the average number of photons in the signal 
is always equal the the mean number of photons in the idler. 

From the time dependent bosonic operators given in eqs.(|5]) and ([6]) we 
compute the average photon number in mode j(= 1, 2), 

(%(^)) = -r^\C\xr) - y^)n^, + S\xT){n,^,,, + 1)] (9) 



Because the pump amplitude A, is treated as constant, the solution to the 
parametric interaction cease to be valid when the number of pair of photons 
is such that determines an appreciable depletion of the pump mode, this 
limits our calculation to the case {fijij)) « |^p and t = gt < 1. 
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Figure 1: (a) Time- dependent average number of photons (uj) in mode j = 
1, 2 for uncorrelated vacuum inputs: y = (thick line), y = 0.9 (dotted line). 
(b) Average number of photons as function of y, for r = 0.9. 

For small interaction times r << 1, C'^{xt) ~ 1, and the number of 
photons produced is insensitive to the mismatch y. For increasing interaction 
times the average number of photons in each mode grows exponentially with 
C'^{xt) and S'^{xt). For a given interaction time the squeezing parameter 
r{y) is smaller for bigger mismatch, and we expect smaller production of 
photons with increasing mismatch. These conclusions are corroborated by 
the plots of Figure 1. 

In Fig.l (a) is plotted the time- dependent average number of photons in 
the signal for nio = ?t.2o = 0: y = (thick line) and y = 0.9 (dotted line). 
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In Fig.l (b) the the average number of photons is plotted as function of the 
mismatch parameter y, for r = 0.9. These graphics show that the average 
number of photons produced by the nonhnear interaction increases with the 
interaction time r, and decreases with increasing mismatch y. 

3 Degree of entanglement 

The concepts of entanglement and information encoded in a quantum state 
are closely related. For pure states, bipartite entanglement, is equivalent to 
the lack of information (mixedness) of the reduced state of each subsystem. 
For mixed states, each subsystem has its own level of impurity, and moreover 
the global state is itself characterized by a nonzero mixedness which makes 
difficult to define a single measure of entanglement. 

We assume that our system can be in three different initial states. In the 
first case, modes 1 and 2 are excited from the vacuum state, i.e., the initial 
state is |0)i|0)2, which is the case of the spontaneous parametric downcon- 
version. In the second case the initial state is the product of a coherent state 
|a)i, of the injected signal, by the vacuum state |0)2, for the idler, i.e., the 
initial state is |a)i|0)2. These two quantum initial states are pure states and 
the unitary evolution of the system will keep them pure states. 

The third initial state that we consider is a mixed state given by the 
product of two thermal states in modes 1 and 2, 

PO = Pio ® p20 (10) 

where 

CXD 

p.o= Ee-^^"^>.)(^.l (11) 

nj=0 

is the density operator of mode j(= 1, 2) in thermal equilibrium at tempera- 
ture T, with Pj = -j^ and ks is the Boltzmann constant. The initial average 
number of photons in mode j is tijq = -j^ — . Though at room temperature 
fijQ << 1 in the optical part of the electromagnetic spectrum, in the mi- 
crowave part of the spectrum fijQ » 1 and we cannot ignore the presence 
of thermal photons in the nonlinear medium. 

All these three initial quantum states are Gaussian states, i.e. they have 
Gaussian Wigner and characteristic functions and hence they are fully char- 
acterized by the first and second statistical moments [131 E] of the dimension- 



^ (cLi—dy 



less mode conjugate quadratures Xj = y ^(% + o,j) and pj 



(j = 1, 2), that is, by the vector of the mean values X = ((xi), (£2), {pi), {P2)) 
and by the covariance matrix (CM) a 



cr. 



V 



{xiXj + XjXi) - {xi){xj) (12) 



The first moments can be arbitrarily adjusted by local unitary operations. 
Such operations leave entropy and entanglement invariant. For uncorrelated 
vacuum inputs, the first moments are already zero, without any need of 
adjustment. 

The correlations between the quadratures of modes 1 and 2, are encoded 
in the off diagonal blocks of the CM 

7=f"" "") (13) 

\ 0-23 0-24/ 

Unitary evolutions, generated by bilinear Hamiltonians such as ([1]), pre- 
serve the Gaussian statistics of the quantum states fTH [T5j , hence the quan- 
tum state of the coupled modes at any instant t > 0, is still a Gaussian 
state. 

The degree of entanglement in a pure bipartite system is equivalent to 
the degree of mixedness of each subsystem and is properly quantified by 
the entropy of entanglement S, defined as the Von Neumann entropy of its 
reduced density operator [131 [H] . In terms of the entries of the CM a, the 
Von Neumann entropy of each mode is given by [12] . 

£ = /(^i) (14) 

and 

/(^) = (^ + 2) 1^(^ + 2) - (^ - 2^ ^""^^ ~ 2^ ^^^^ 



where Uj = Ja-jjaj+ij+i — cr'j^j+i > 1/2 is the symplectic eigenvalue of the 
reduced single- mode CM o"j (j = 1,2). Each mode displays maximal disorder 
subjected to given physical constrains. Our system is unbounded, therefore 
the entropy is constrained by specifying its mean energy. Whenever Vj = 1/2, 
mode j is in a pure state and the composite system is in a separable quantum 
state. 



For two uncorrelated vacuum inputs |0)i|0)2, the initial average values of 
photons in each mode is uiq = 77,20 = and for the initial state |a)i|0)2 they 
are tt-io = |ap and 77,20 = and the symplectic eigenvalues of the reduced 
covariance matrices cxj are given by 




^^1 = ^^2 = 7T^/ 1 + 4 ; .2 i_ 2 (2^10 + 1) (16) 



I C'^{xt) — 7/2 S'^{xt) 

1-7/2 1-7/2 

showing that Uj > 1/2 for t > 0, then the entropy of each mode is greater 
than zero and the two modes are entangled from the beginning of the interac- 
tion. For interaction times such that the squeezing parameter is r{y) << 1, 
the entropy of each mode is almost insensitive to the mismatch, for increasing 
interaction times the exponential growth of the hyperbolic functions domi- 
nates and the degree of entanglement decreases with increasing mismatch. 
This behavior is expressed in Figure 2. 

Figure 2(a) shows the time-dependent degree of entanglement S for the 
initial state |0)i|0)2, for zero mismatch (7/ = 0) (thick line) and for 7/ = 0.9 
(dotted line). In panel (b) the degree of entanglement is plotted as function 
of the dimensionless mismatch 7/ for a given dimensionless interaction time, 
r = 0.9. As we expected, the degree of entanglement increases with the in- 
teraction time and, for a given interaction time decreases when the mismatch 
increases. 

In mixed states the mere distinction between classical correlations, i.e. 
producible by means of local operations and classical communication (LOCC) 
only, and entanglement, due to a purely quantum interaction between sub- 
systems, is a highly nontrivial task and the entropies of each mode are no 
more a convenient measure of entanglement. 

In our bipartite system, with two interacting modes, we can use the partial 
transpose (PPT) criterion which says that a quantum state is separable if 
and only if i>_ > 1/2, where i>_ is the smallest symplectic eigenvalue of the 
partial transpose of a, given by [I2] 




i>^ = J-{A{cx)-JA{a)-4deta) (17) 



where A(cr) = det ai + det (72 — 2 det 7. The logarithmic negativity provides 
a proper quantification of entanglement in two-mode Gaussian states and is 
a computable measure of entanglement, defined by [161 [HI 

^ _ f -ln2i>_ if -ln2z>_>0 .^„. 

^^-\ if -ln2z>_<0 (^^^ 
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Figure 2: (a) Time- dependent degree of entanglement £ for the initial state 
|0)i|0)2 : y = (thick line) and y = 0.9 (dotted line), (b) S is plotted as 
function of y for r = 0.9. 

For the initial uncorrelated thermal state given in f IlOp the entries of the 
covariant matrix are given in the Appendix and are used to compute the 
logarithmic negativity Sj\f which is plotted in Figure 3. 

Figure 3 (a) shows, Sj^, as function of time: dashed line corresponds to 
uncorrelated vacuum input state with y = and dotted line is the same input 
state at y = 0.9, thick line corresponds to the initial thermal state for y = 
and the dasheddot line is the same initial condition for y = 0.9. The function 
£^ is always positive for uncorrelated vacuum inputs, independently of the 
mismatch parameter and for any interaction time, in accordance with our 
previous conclusions. For the initial thermal state, Sj\f becomes positive only 
after a certain time of interaction. This clearly shows that the entanglement 
is more difficult to attain for a mixed thermal state. Moreover, the com- 
posite system needs more time of interaction to become entangled when the 
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mismatch is nonzero. 

In Fig. 3 (b) we plot £j^ for the case of the initial thermal state, as function 
of the mismatch and for r = 0.7. For this interaction time, £jsf is always 
positive, decreasing with the mismatch meaning that the system is entangled. 

Finally, in Fig. 3(c) we plot different quantities for the initial thermal 
state with tt-iq = 1, ^20 = 2 : the Von Neumann entropy of mode 1 for y = 
(dashed line) and for y = 0.9 (dotted line), and Sj^ for y = 0.9 (thick line). 
The entropy of mode 1 increases with time and, for zero mismatch, it is 
always greater than for nonzero mismatch. 




Figure 3: (a) Sj\f as function of time for the uncorrelated vacuum input state: 
y = (dashed line) and y = 0.9 (dotted line). The same quantity for the 
initial thermal: y = (thick line) and y = 0.9 (dot-dashed line) is the same 
initial condition for y = 0.9. (b) Sj\f as function of y for r = 0.7. 

The degree of entanglement for all three initial states that we have stud- 
ied, decreases with increasing mismatch. This is a direct consequence of the 
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Figure 4: (c) Von Neumann entropy of mode 1: y = (dashed line) and 
y = 0.9 (dotted line), and Sj^ for y = 0.9 (thick line). All three functions 
were computed for the initial thermal state [ld\) with nio = l,n2o = 2. 



dependence of the effective squeezing parameter r{y), present in the time 
dependent bosonic operators ([5]) and (Q, with the phase-mismatch. 

4 Summary and concluding remarks 

We have extended the computation of the degree of entanglement, in non- 
degerate two-mode parametric processes in order to include the effect of a 
phase mismatch. A correction factor was introduced in the usual hamil- 
tonean of these processes and exact solutions of the Heisenberg equations 
were derived. 

It was shown that the new squeezing parameter includes a dependence 
in the mismatch parameter and that its intensity decreases with increasing 
mismatch. This behavior is responsible for generating an average number of 
photons in both modes that decreases with increasing mismatch. 

The degree of entanglement was computed for three different initial quan- 
tum states. The two pure initial states are associated with the downconver- 
sion and with the injected signal of the amplifier. The initial mixed thermal 
sate takes into account the finite temperature of the nonlinear medium and 
may be relevant for processes in the microwave region of the electromagnetic 
spectrum. The degree of entanglement shows a similar behavior in these three 
cases, i.e, for a given interaction time it decreases with increasing mismatch 
and increases with the generated number of photon pairs. 
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To measure the degree of entanglement that we have derived in this paper 
we propose to make measurements of the statistical properties of the quadra- 
tures correlation of the two-modes, which are the entries of the covariance 
matrix of the system, trough the usual techniques of homodyne detection 
[I81[I9]. 

The Hamiltonean including the phase-mismatch factor is similar to the 
Hamiltonean associated with a slight detuning between the pump, the signal 
and the idler, in other words this Hamiltonean takes into account a possi- 
ble effect of the finite bandwidth of the modes in interaction. In this case, 
our results suggest a degradation of the quality of entanglement when finite 
bandwidth of the modes cannot be ignored. 



5 Appendix 

Using the time dependent bosonic operators given in eqs.([5]) and ([6]) we 
compute the entries of the covariance matrix a for the thermal initial state 
(TTUj) . obtaining 



2 

0-12 (t) = 0-21 (r) = 0-34 (r) = 0-43 (r) = 

0-24 (r) = 0-13 (r) ; 0-23 (r) = 







(J = l,2) 



0-14 (r) = (nio + ri2o + 1 



cri4lT-j 



y 



''^^Ci.r)Si.r) 



(19) 

(20) 
(21) 

(22) 



sm{2uj'xT)y 2 



y2 



S'ixr] 



(23) 
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Figures captions 

Fig. 2 (a) Time-dependent degree of entanglement £ for tlie initial state 
|0)i|0)2 : y = (thick line) and y = 0.9 (dotted line), (b) £ is plotted 
as function of y for r = 0.9. 

Fig. 3 (a) £f^ as function of time for the uncorrelated vacuum input 
state: y = (dashed line) and y = 0.9 (dotted line). The same quantity 
for the initial thermal: y = (thick line) and y = 0.9 (dot-dashed line) 
is the same initial condition for y = 0.9. (b) £j\/ as function of y for 
r = 0.7. (c) Von Neumann entropy of mode 1: y = (dashed line) and 
y = 0.9 (dotted line), and £^f for y = 0.9 (thick line). All three functions 
were computed for the initial thermal state (fTUj) with uiq = 1, ^20 = 2. 
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